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Abstract. This letter treats the quantum random walk on the line determined by a 2 x 2 
unitary matrix U. A combinatorial expression for the mth moment of the quantum random walk is 
presented by using 4 matrices, P, Q, R and S given by U. The dependence of the mth moment on 
U and initial qubit state <p is clarified. A new type of limit theorems for the quantum walk is given. 
Furthermore necessary and sufficient conditions for symmetry of distribution for the quantum walk 
is presented. Our results show that the behavior of quantum random walk is striking different from 
that of the classical ramdom walk. 



1 Introduction 

Recently various problems for quantum random walks have been widely investigated by a 
number of groups in connection with quantum computing. Examples include Aharonov 
et al.^ 1 ' Ambainis et ai, ^ Bach et al.,^ Childs, Farhi and Gutmann/ 4 ) Diir et al.,( 5 ' 
Kempe/ 6 -*, Konno/ 7 ) Konno, Namiki and Soshi/ 8 ) Konno et a/./ 9 -* Mackay et a/./ 10 ) Moore 
and Russell/ 11 ^ Travaglione and Milburn/ 12 ) Yamasaki, Kobayashi and Imai.*- 13 - 1 For a more 
general setting including quantum cellular automata, see Meyer. ( 14 ) A more mathematical 
point of view of quantum computing can be found in Brylinsky and Client 15 ) 

In Ambainis et al.,^ they gave two general ideas for analyzing quantum random walks. 
One is the path integral approach, the other is the Schrodinger approach. In this paper, we 
take the path integral approach, that is, the probability amplitude of a state for the quantum 
random walk is given as a combinatorial sum over all possible paths leading to that state. 

The quantum random walk considered here is determined by a 2 x 2 unitary matrix U 
stated below. The new points of this paper is to introduce 4 matrices, P,Q,R, S given by 
the unitary matrix U, to obtain a combinatorial expression for the characteristic function by 
using them, and to clarify a dependence of the mth moment of the quantum random walk on 
the unitary matrix U and initial qubit state ip. We give a new type of limit theorems for by 
using our results. Moreover we present necessary and sufficient conditions for symmetry of 
distribution for the quantum walk. Our results show that the behavior of quantum random 
walk is remarkable different from that of the classical ramdom walk. 
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The time evolution of the one- dimensional quantum random walk studied here is given 
by the following unitary matrix (see Nielsen and Chuang^ 16 )): 



U = 



a b 
c d 



c| 2 + 



where a,b,c,d G C and C is the set of complex numbers. So we have \a\ 2 + \b\ 2 
\d\ 2 — 1, ac + bd — 0, c = — Ab, d = Aa, where z is a complex conjugate of z G C and 
A = det U = ad — be with |A| = 1. The qunatum walk is a quantum generalization of 
the classical random walk with an additional degree of freedom called the chirality. The 
chirality takes values left and right, and means the direction of the motion of the particle. 
The evolution of the qunatum walk is given by the following way. At each time step, if the 
particle has the left chirality, it moves one step to the left, and if it has the right chirality, it 
moves one step to the right. The unitary matrix U acts on two chirality states \L) and \R): 



\R) 



a\L) +c\R) 
b\L)+d\R) 



where L and R refer to the right and left chirality state respectively. In fact, define 



1 




\R) 




1 



so we have 



U\L) = a\L) + c\R) 
U\R) = b\L)+d\R) 

More precisely, at any time n, the amplitude of the location of the particle is defined by 
a 2-vector G C 2 at each location {. . . , —2, —1, 0, 1, 2, . . .}. The probability the particle is at 
location k is given by the square of the modulus of the vector at k. If (^(n)) define the 
amplitude at time n at location k where 



l*kW> = 



with the chirality being left (upper component) or right (lower component), then the dy- 
namics for \^ k (n)) in quantum walks is given by the following transformation: 

|tt fc (n + 1)) = \L){L\U\* k+1 (n)) + \R)(R\U\^ k ^(n)) 

The above equation can be rewritten as 

|* fc (n + 1)) = P|* fc+ i(n)> + Q|tt fc -i(n)) 

where 



P = 



a b 




Q 





c d 
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We see that U = P + Q. The unitarity of U ensures that the amplitude always defines a 
probability distribution for the location. 

The simplest and well-studied example of a qunatum walk is the Hadamard walk whose 
unitary matrix U is defined by 



H 



1 

72 



1 1 
1 -1 



The dynamics of this walk corresponds to that of the symmetric random walk in the classical 
case. In general, the following unitary matrices can also lead to symmetric walks: 



IT] 



77,0,1/1 



e 
V2 



-i(4>+ip) 



where 77, 0, and ip are real numbers (see, for example, pp. 175-176 in Nielsen and Chuang^ 16 )). 
In particular, we see C/0,0,0 = H. However symmetry of the Hadamard walk depends heavily 
on the initial qubit state, see Theorem 4 in this letter or Ref. 8. 

In the present letter, the study on the dependence of some important quantities (e.g., 
the mth moment, limit distribution) on initial qubit state is one of the essential parts, so we 
define the set of initial qubit states as follows: 



GC 2 :|a| 2 + |/3| 2 = l 





a 







The next section is devoted to a method given by us and some our new results. 



2 Our Method and Results 

Let X% be probability distribution of the quantum random walk starting from initial qubit 
state if e $. Note that P(X$ = 0) = 1. In order to study P(X% = k) for n + k = even, we 
need to know the following expression. For fixed I and m with I + m = n and m — I = k, 

E(l, m)=J2 P h Q mi P h Q m2 ■ ■ ■ P ln Q mn 

summed over all lj, rrij > satisfying mi + • • • + m n = m and l\ + • • • + l n — I. Note that 

P(X* = k) = m,m)<pYm,m)<p) 

where * means the adjoint operator. For example, in the case of P(X% = —2), we have to 
know the expression, 5(3, 1) = P 3 Q + P 2 QP + PQP 2 + QP 3 . By using P 2 = aP, Q 2 = dQ, 
the above equation becomes S(3, 1) = a 2 PQ + aPQP + aPQP + d 2 QP. In our treatment of 
quantum walks, as well as the matrices P and Q, it is convenient to introduce 



R = 



c d 




S 





a b 



The next table of products of P, Q, R, and S is very useful in computing some quantities: 
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P 


Q 


R 


S 


p 


aP 


bR 


aR 


bP 


Q 


cS 


dQ 


cQ 


dS 


R 


cP 


dR 


cR 


dP 


S 


aS 


bQ 


aQ 


bS 



where PQ = bR. We should remark that P,Q,R, and S form an orthonormal basis of 
the vector space of complex 2x2 matrices M 2 (C) with respect to the trace inner product 
(A\B) = tr (A*B). So Z(l,m) has the following form: 

E(l, m) = p n (l, m)P + q n (l, m)Q + r n (l, m)R + s n (l, m)S 

Next problem is to obtain explicit forms of p n (l,m), q n (l,m),r n (l,m) and s n (l,m). For 
example, we have S(3, 1) = labcP + a 2 bR + a 2 cS, therefore, Pa(3, 1) = 2abc, 54(3, 1) = 
0, r 4 (3, 1) = a 2 b, 54(3, 1) = a 2 c. In abed = case, the argument is much easier. So from 
now on we focus only on abed 7^ case. In this case, the next key lemma is obtained by a 
combinatorial method. 



LEMMA 1 We consider quantum random walks in one dimension with abed 7^ 0. Suppose 
that I, m > with I + m = n, then we have 

(i) for I A m(= min{/, m}) > 1, 

(ii) for /(= n) > 1, m — ; 
(Hi) for I = 0, m{= n) > 1, 

By this lemma, the characteristic function of X% for abed 7^ case is obtained. Moreover, 
the mth moment of X% can be also derived from the characteristic function in the standard 
fashion. Here we give only the result of the mth moment. Result on the characteristic 
function and proofs of our results mentioned here appear in our mathematical paper/ 7 ) 



07 Aery Ao b 



E(l,0) = a l - l P 
5(0, m) = A^a^Q 



PROPOSITION 2 We consider quantum random walks with abed 7^ 0. 
(i) When m is odd, we have 



e((x:d = \ a 



2(71-1) 



-n 



[^2~\ k k 

+ E EE 

k=l 7=15=1 



a\ 2 - \b\ 2 ) (\a\ 2 - \{3\ 2 ) +2(aab[3 + aab[3)} 



\b\ 2 Y +S f k — l\ f k — l\ fn — k — l\ fn — k — V 



v 7 — 1 J \5 — ly v 7 — 1 



5-1 



X 



(n - 2k) m+1 

•y5 



-{n(\a\ 2 -\b\ 2 )+ 1 + 5}(\a\ 2 -\{3\ 2 ) 
1 + 5 



+ 



2n (aab(3 + aab(3) 



(ii) When m is even, we have 
E((X%) m ) = 



n 



V2 \ k k 

+ E EE 

k=l 7=1,5=1 

in - 2k) 1 



\b\ 2 V +S fk - 1\ fk - l\ fn - k - l\ fn - k - T 



x 



7<5 



n 



\a\ 2 J \7-lJ\5-lJ\ 
-k) 2 + k 2 -n{ 1 + 5) + ^ 



5-1 



It should be noted that for any case, when m is even, E{{X^) m ) is independent of initial 
qubit state ip. Moreover we have the following new type of limit theorems: 



THEOREM 3 We assume abed ^ 0. If n -> 00 , then 



n 

where Z v has a density 



for x G (— \a\, \a\) with 

E(Zn = - ( \a\ 2 - W + + x (i _ ^H 2 ) 



e{{z*) 2 ) = 1 - v/i-M 2 

and Y n =>- Y means that Y n converges in distribution to a limit Y . 



We remark that standard deviation of Z v is not independent of initial qubit state p = 
/?]. The above limit theorem suggests the following result on symmetry of distribution 
for quantum random walks. This is a generalization of the result give by Konno, Namiki 
and Soshi® for the Hadamard walk. Define 

$ s = {p g $ : P(X% = k) = P(X% = -k) for any n G Z+ and k G Z} 

$0 = {(/? G $ : E(X%) = for any n G Z+} 

$_l = {</? = G $ : \a\ = \[3\, aabfl + aabf3 = oj 
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and Z (resp. Z + ) is the set of (resp. non-negative) integers. For (p G $ s , the probability 
distribution of X% is symmetric for any n G Z + . Using the explicit form of E(X%) given by 
Proposition 2 (i) (m = 1 case), we have 



THEOREM 4 We assume abed ^ 0. Then we have 

$ s = $ = $j. 

From now on we give an outline of our proof of Theorem 3 (for more details, see Konno^). 
To do so, we introduce the Jacobi polynomial P^(x), where P^(x) is orthogonal on [—1, 1] 
with respect to (1 — x) u (l + xY with v, /j, > —1. 

First we obtain the next asymptotics of characteristic function E(e l ^ Xn ^ n ) by using the 
Jacobi polynomial: if n — > oo with k/n = x G (—(1 — M)/2, (1 + \a\)/2), then 

E{e* x *' n ) ~ £ |a| 2 "- 4fc ->| 4 
fc=i 



x 



2 " 2 J a+1 (^r afc ) 2 - f^r 2 *^- 5 * + ^(n°iv 2fe ) 2 } cos((i - 2x)o 

! " " r --{(|a| 2 - \b\ 2 ){\a\ 2 - \P\ 2 ) + 2{aa^ + aab(3)}{P^ 1 - 2k ) 2 



x /IX 



oil |2 |o|2 + aab[3 \ D 0,n-2fc D l,n-2fcl • • //-, Nt\ 

-2 Mai j^p ) P fc-i P fc-i }«sm((l - 2x)0 

where /(n) ~ g(n) means f(n)/g(n) — > l(n — > oo), and P^^ 2fe = P^^ 2fc (2|a| 2 — 1) (i = 0, 1). 

Next we prepare the following asymptotic results for the Jacobi polynomial p®+ an fi+ bn ( x ) 
derived by Chen and Ismail/ 17 ) : if n — > oo with fc/ra = rr G (—(1 — H)/2, (1 + \a\)/2), then 



fell ~ - - 



2k-n 

cos(An + 5) 



|2fc-n 



7m 



o _ n / ~ 

~7fe ^-^-w - (/tn+B+fl) 

where A = (1 — |a| 2 )(4x 2 — 4x + 1 — |a| 2 ), A and -B are some constants (which are independent 
of n), and 9 G [0, 7r/2] is determined by cos 9 = ^(l-\a\ 2 )/Ax(l-x). 
Combining these results, Theorem 3 is obtained. 

Now we compare our analytical results (Theorem 3) with the numerical ones* 10 ' 12 -' for the 
the Haramard walk. In this case, Theorem 3 implies that if — v / 2/2 < a < b < a/2/2, then 
as n — > oo, 



P(a < XII n < b) -> f 

J a 



1 - (H 2 - \(3\ 2 + a(3 + a(3)x , 

^ dec 

7i(l-x 2 )Vl-2x 2 
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for any initial qubit state tp = t [a,f3}. For the classical symmetric random walk Y° starting 
from the origin, the well-known central limit theorem implies that if — oo < a < b < oo, then 

as n — > oo, 




This result is often called the de Moivre-Laplace theorem. When we take y? = ^1 /y/2, if y/2] 
(symmetric case), then we have the following quantum version of the de Moivre-Laplace 
theorem: if — \[2j2 < a < b < y/2/2, then as n — > oo, 

P(a < Xtln <b)^ [ b 1 - dx 

So there is a remarkable difference between the quantum random walk X% and the classical 
one Y° even in a symmetric case for ip = t [l/y/2,i/y/2]. 
Noting that E{X£) = (n > 0) for any ip G we have 

sd(X%)/n -> j(2-\/2)/2 = 0.54119 . . . 

where sd(X) is the standard deviation of X. This rigorous result reveals that numerical 
simulation result 3/5 = 0.6 given by Travaglione and Milburn^ 12 ) is not so accurate. 

As in a similar way, when we take (p = t [0,e td ] where 6 e [0, 2tt) (asymmetric case), we 
see that if —y/2/2 < a < b < y/2/2, then as n — > oo, 

P(a < XHn < b) - f — ly—^ dx 

Ja 7r(l — x)vl — 2a; 2 

So we have 

E(X%)/n -> (2 - v^)/2 = 0.29289 . . . , sd(X%)/n -> j(y/2-l)/2 = 0.45508 . . . 

When = *[0, 1] (6» = 0), Ambainis e* a/. (2) gave the same result. In the paper, they took 
two approaches, that is, the Schrodinger approach and the path integral approach. However 
their result comes mainly from the Schrodinger approach by using a Fourier analysis. The 
details on the derivation based on the path integral approach is not so clear compared with 
Konno( 7 ). 

In another asymmetric case ip = t [e te ,0] where G [0, 2ir), a similar argument implies 
that if -y/2/2 < a <b < V2/2, then as n -> oo, 

P{a < Xljn < b) - f 1 dx 

Ja 7T(1 + XjV 1 — 2x l 

Note that /(— x; 1 [e ie , 0]) = /(x; *[0, e* e ]) for any x G (-^2/2,^2/2). Therefore concerning 
the mth moment of the limit distribution, we have the same result as in the previous case ip = 

t [0,e te ]. So the standard deviation of the limit distribution Z v is given by \J (a/2 — l)/2 = 
0.45508 .... Simulation result 0.4544 ± 0.0012 in Mackay et alS 10) (their case is 9 = 0) is 
consistent with our rigorous result. 
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